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Abstract 

We present a model of bi-phasic vesicle in the limit of large surface 
tension. In this regime, the vesicle is completely stretched and well 
described by two spherical caps with a fold which concentrates the 
membrane stress. The conservation laws and geometric constraints 
restrict the space of possible shapes to a pair of solutions labeled by a 
parameter r given by line tension/pressure. For a given r value, the 
two solutions differ by the length of the interface between domains. 
For a critical value r c , the two vesicle shapes become identical and no 
solution exists above this critical value. This model sheds new light on 
two proposed mechanisms (osmotic shocks and molecule absorption) 
to explain the budding and the fission in recent experiments. 
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1 Introduction 



The cell membrane is a bilayer made out of a mixture of lipid species. The 
membrane is both the boundary of the cell and an interface inside the cell, 
separating different compartments. This soft structure is responsible for 
many biological properties. Intracellular traffic is also realized by mem- 
brane structures: a membrane vesicle buds from one compartment, trav- 
els through the cytosol and fuses with another compartment. Despite the 
fluidity of the lipid bilayer, the cellular membrane presents a lateral in- 
homogeneity due the formation of dynamical microdomains, called rafts 
( Simons and Ikonen, 1997| ). These microdomains have been shown to be rich 



in cholesterol and sphingolipid (B rown and London, 2000| ). In vivo, the rafts 
have not been directly observed but their size has been estimated to be be- 
tween 20 and 700 nm dChazal and Gerlier, 2003 1 . A central question in mem- 
brane biology and biophysics is to understand how this spatial organization 
is used by the cell, in particular to favor interactions with proteins. Due to 
their size and specific composition, it has been argued that rafts play a role in 



protein docking, signaling, intracellular traffic (van Meer and Sprong, 2004) 



or virus budding dChazal and Gerlier, 2 003). 

Recently, a model system of Giant Unilammelar Vesicles (GUV) includ- 
ing sphingomyelin-cholesterol domains was developed dDietrich et al., 2001 ). 
These domains, which are supposed to reproduce raft composition, are the re- 
sult of a phase separation of the lipid species dVeatch and Keller, 2003D . They 
are more structured than the surrounding classical liquid bilayer but remain 
in a liquid state. For this reason, they are called "liquid-ordered" domains 
whereas the classical membrane is called "liquid-disordered". A large num- 
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ber of studies have focused on the thermodynamic of liquid-ordered phases, 
in particular the effect of temperature or composition changes on domain 



formation flde Almeida et al., 2003} |Veatch an d Keller, 2003). Multi-phase 
vesicles are elegant and efficient tools to study the mechanical properties of 
microdomains. It can be used to understand how rafts bud and make daugh- 
ter vesicles for intracellular traffic, but also how detergent addition can isolate 
rafts from the cell membrane. Recent experiments have shown that liquid- 
ordered domains can be separated from the initial vesicle by using tubu- 



lar deformations JAilain et al., 2004| ), osmotic shocks (Baum gart et al., 2 003 



Bassereau and Roux, personnal communication! or absorption of external molecules 



like proteins or detergents ( |Staneva et al., 2004 |Staneva et al., submitted ). 
Here, we develop a macroscopic theory for the two last situations. Our model 
describes the budding preceeding the fission where the liquid-ordered domain 
is lift up from the liquid-disordered vesicle. 

Budding and fission have already attracted many theoretical works for ho- 
mogeneous ( jJaric et al., 1995[[Seifert, 1997t|Do"bereiner et al., 1997[|Tanaka et al., 2004] 
Sens, 2004) or inhomogeneous ( Seifert, 1993 Jtilicher and Lipowsky, 1 996 



|Kohyama et al., 2003HLaradji and Sunil Kumar, 2004 Harden et al., submitted 1 
membranes. The models vary depending on the physical interactions in- 
volved but they are all based on the minimization of the bilayer energy 
( Helfrich, 1973D - Due to the non-linearity of the steady-shape equations, a 
numerical treatment is often required. We focus our attention on multi-phase 
vesicles slightly stretched, a situation often encountered in experiments. In 
this case, osmotic pressure effects dominate and we show that the vesicle can 
be described by two spherical caps with an elastic junction. The variational 
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procedure with surface constraints allows to find two solutions for any ratio 
t = line tension/pressure less than a critical value r c . The stable solution is 
the one observed experimentally. An osmotic shock increases the control pa- 
rameter t and so destabilizes the stable solution which may lead the system 
to a complete fission of the neck. The case of detergents is slightly different 
since it requires an energy model for molecular absorption in the membrane. 
When detergent molecules are added in the membrane, they locally deform 
the bilayer. According to Leibler's model ( Leibler, 1986| ), these curvature 
defects can be taken into account by a term in the energy proportional to 
both the average curvature and the concentration of molecules. Homoge- 
neous concentration of molecules is favored away from the interface between 
domains. At the junction, a concentration gradient appears. If the chemical 
inhomogeneity is localized at the junction, the addition of molecules leads 
to an increase of the effective line tension inducing a budding and a possible 
separation into two independent vesicles. 

Our model explains qualitatively and even quantitatively the budding 
and fission created by osmotic shocks or proteins absorption. It is a physical 
approach based on continuum description and its domain of validity ends at 
the molecular level. Because of its simplicity, extension and application to 
other processes may be achieved easily. 
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2 Membrane description. 



2.1 Inhomogeneous lipid bilayer. 

We consider an inhomogeneous vesicle constituted by two lipid phases: a 
'liquid-ordered' (l ) and a 'liquid-disordered' (la). Both phases are in the 
liquid state but the (l ) domain is more structured than the (Id) phase due 
to the following reasons: there are specific interactions between molecules 



( Li et al., 2001 ) and/or there is an optimization of biphilic space packing 
dHolopainen et al., 2004D . Steady morphologies and their out-of-plane defor- 
mations are well described by the Canham and Helfrich's model with energy 



for each phase i given by: 



2/t i H 2 + ^ ) K + S i 



dS 



s 



H and K are respectively the mean and Gaussian curvature. The elastic bend- 
ing rigidity «j and Gaussian rigidity Kq are expected to be higher in the l Q 



phase. Typical values can be found for example in (Lipows ky and Sackmann, 1995 ) 
Ki d ~ 20kbT and Ki o ~ 80kbT. Values of Gaussian moduli are notori- 
ously more difficult to measure but a recent study mentions values of order 



Kq = — 0.83Kj dSiegel and Kozlov, 2004| ). Although is a surface integral 



the Gaussian contribution to the energy is indeed a contour integral calcu- 
lated at the interface between the two domains, due to the Gauss-Bonnet 
theorem. 

The last contribution in EqJT] is related to the possible extension of the 
membrane. In the case of a stretched vesicle, this contribution is large com- 
pared to the elastic energy and the membrane surface can be considered as 
constant. This is taken into account by introducing the Lagrange multiplier 



The total energy of the two-domain vesicle includes the energy JH) of each 
phase plus two coupling terms. First, a sharp interface of vanishing thickness 
exists between the l and Id phases. Any increase of its length requires an 
energy proportional to a line tension a. Second, the vesicle membrane is 
lightly permeable to the water but not to the ions or big molecules present 
in the surrounding water medium. This induces an osmotic pressure P. The 
energy of the coupling terms is: 



c 

2.2 Proteins or detergent-membrane interactions 

External molecules such as proteins or detergents can be absorbed in both 
phases but with different efficiencies. Their introduction in the membrane is 
well described by a Landau's approach with an optimal homogeneous con- 
centration <ft eq . Departure from this value or inhomogeneity of concentration 
<fi has a cost in energy, assumed quadratic to leading order. The energy cost 
is given by two positive constants in each phase: and If the proteins are 
soluble or not in the surrounding medium, we can either set the number of 
these molecules in each phase or set the chemical potential /ij of the phase i. 
We choose to fix /ij but this has no real incidence on the results since it only 
affects the definition of /ij. Therefore the free chemical energy of absorption 
for each phase is: 



f p = / (4(<l> - + ^ W) 2 + ^)ds+ amhs (3) 




(2) 
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The last integral in EqOl represents the local distortion of the membrane 



induced by the absorbed molecules (Leibler, 1986 Bickel et al., 2001 ). It is 



proportional to the mean curvature of the membrane with a weight depending 
on the local concentration 0, as suggested by S. Leibler (Leibler, 1986), Aj 



being a coupling constant. The absorption process itself affects differently 
the two leaflets of the vesicle. We restrict our attention to the case where 
the adsorption takes place on one side only. In such case, Aj is positive on 
the outer side absorption and negative on the inner side. When the two 
layers are affected by absorption, two concentration fields are necessary and 
our theoretical framework can be easily adapted to address such situation. 
Taken into account all previous contributions, the total free energy for the 
system is given by: 

F TOT = F° m + F ^ + F; + F p d + F c (4) 

The usual variation procedure to identify extrema of this energy produces 
the so called Euler-Lagrange equations. 



2.3 Euler-Lagrange equations. 

Minimization of the free energy Ftot gives the static solutions for the mem- 
brane. Looking for axisymmetric shapes, we choose the cylindrical coordi- 
nates and we parameterize the surface by the arc-length s. The vesicle shape 
is given by r(s) and i/j(s) (see Fig. [TJ. We have derived the Euler-Lagrange 
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equations associated with F TOT ( Allain and Ben Amar, 2004| ). They are 



ship) cos(^) f Pr 7 . A; 
5 cos(^) - - — coshp) H sin ^ H 



i = ji/> ,1 -—^) 2 + t i -Pr^)-k l ^' + ^<p 2 (5b) 

a" -/CosW Ai /sin(^) A on hi 

r' = cos(ip). (5d) 

These equations have to be solved with the suitable boundary conditions at 
the border between the two domains. To simplify the notations, we introduce 
the following parameters: = Sj + a i /20g ? . and jli = ^i — ai4> eqi . Assuming 
continuity of both the radius r, the angle ip and the molecules concentra- 
tion 0, the variational procedure gives also three equations for the boundary 
conditions: 

Kiip'{sj - e)r(sj) + («i + K Gl ) sin(V>(s j)) - A^i-s j)r(sj) (6a) 

-K 2 ip'(sj + e)r(sj) - (k 2 + k G2 ) sin(^(sj)) + A 2 <fr(sj)r(sj) = 0, 
l(s J -e)- 1 (sj + e) + a = (6b) 
M(8j - e) - (3 2 <P\sj + e) = 0. (6c) 

where sj is the arc-length at the junction, label 1 denotes the phase for 
s < sj and label 2 the phase for s > sj. 

Since these equations are highly non-linear, there is no exact solutions 
for arbitrary values of the physical parameters. However further analyti- 
cal progress can be obtained in the limit of large pressure (stretched mem- 
brane). Remarkably, this treatment only requires simple analytical alge- 
bra and allows to explain experimental features such as the complete bud- 
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ding of the ordered phase obtained by different groups using either osmotic 
shocks (Baumg art et al., 20 03), proteins ( Staneva et al., 2004 1 or detergent 
molecules (Staneva et al., submitted) 



3 Analytical treatment of the membrane shape. 

We first consider a membrane without absorbed molecules. A solution of the 
Euler-Lagrange equations can be easily found if we discard the contribution 
from the elasticity. We use this simple solution as zeroth order and correct it 
for small but not vanishing values of the bending rigidity by using boundary 
layer analysis. We consider also the inclusion of molecules with no chemical 
activity. They are described in the model by curvature defects. For a weak 
coupling between curvature and concentration, a similar strategy is used to 
understand how the molecules affect the membrane shape. 

3.1 The exact zero-order model: the capillary solution. 

For stretched membrane without absorbed molecules (0 = 0), it is believed 
that after electro-formation of GUV vesicles the osmotic pressure dominates 
the elastic energy. When Ki = in both phases, a solution of the Euler- 
Lagrange equation is made of two spherical caps defined by a set of four 
geometrical parameters: the radii of the two caps Ri, R2 and the two angles 
at the boundary 9\ and O2 (see Fig. EJ). The contact between the two caps 
gives a first continuity relation 

#1 sin 0i = # 2 sin 2 - (7) 
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The Euler-Lagrange equations (EqO|) give the values of the two Lagrange 
multipliers Sj and ji, without direct information on the vesicle shape: 

2S, = PR h (8a) 
PR 2 

7i(s) = sin ip cos (8b) 

the angle ^ being proportional to the arc-length s. Only the boundary con- 
dition and the conservation of the area of each phase give the possibility to 
fix completely the ideal shape. From EqGH we deduce: 

R\ sin #i cos 6i = R\ sin 2 cos 2 — — . (9) 

The shape is controlled by the reduced line tension t = cr/P (homogeneous to 
a surface), which can be adjusted by changing the osmotic pressure. As an ex- 
ample, from the figure (lb) in Baumgart et al.'s work (fBaumg art et al ., 2003), 
reproduced here in figurelU we calculate r = 20.5/xm 2 , from estimated values 
of Ri, R 2 , 0i and 2 . Notice that in Baumgart's work, label 1 correspond to 
the Id domain and label 2 to the l D domain. 

Solving EqEl and for the above r-value gives two possible solutions: 
Ri = 5.30/xm, R 2 = 10.5/xm, 0\ = 1.34 and 2 = 0.514, the measured 
values (Fig. 0(a)) but also R\ = 3.97fim, R 2 = 10.3/xm, 0\ = 1.96 and 
2 = 0.364 for the second solution (Fig. 0(b)). In order to explain why the 
first solution is preferred in the experiment, we calculate the energy which is 
restricted here to two contributions: Ftot = —PV + al with I the perimeter 
of the interface. Using a typical length scale L r = 10/im, it is possible to 
construct the dimensionless energy Ftot = F TO t/{.^PLI). Notice that the 
value of L r does not affect the physics of the problem, it is used only to have 
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dimensionless lengths close to 1. So we obtain 

Ftot = -V/nLl + 2tR 1 sin a /i£ (10) 

which gives respectively (-1.380) compared to (-1.377), and shows that the 
experimental observed solution is stable while the other one is unstable as 
expected. 

A systematic study of the pair of solutions for arbitrary values of r is 
straightforward and the results are presented in figure [HI Figure is a 
classical bifurcation diagram when a pair of solutions appears with opposite 
stability. In this problem, r is the control parameter and the energy Ftot is 
the order parameter. These two solutions differ geometrically, the unstable 
solution presenting a smaller neck compared to the stable one (obvious if 
r = 0) (see FigEb)- As r increases, the two solutions become geometrically 
closer up to a finite value of r c = RiR 2 /2. Above the critical value (r > 
r c ), there is no connected solution of the Euler-Lagrange equations but the 
solution with two separated spheres remains. 

This bifurcation diagram describing change in the topology of budding 
spheres is similar to the one found in the catenoid problem where a soap film 
is fixed on two parallel rings, separated by a small distance d compared to 
the radius of the ring. Two different minimal surfaces (with similar catenoid 
shapes) satisfy the variational equations derived from the capillary energy. 
The difference between these two shapes can be measured by the perimeter 
at mid-distance between the two rings. The catenoid with the smaller neck 
is unstable since its area is larger and, experimentally, the other catenoid is 
observed. However, as it is well known, as the distance d increases, the neck 
size decreases, the catenoid is destroyed and is replaced by two independent 
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disks ( |Ben Amar et al., 1998 1 with topology changes. This geometrical in- 
stability is not reversible. At the fission, the neck of the catenoid is not zero 
but the analytical calculation shows that the two catenoids, the stable and 
the unstable, have the same shape. 

In our case, we are faced with the same type of capillary instability where 
there exist two similar solutions whose stabilities are governed by the energy. 
As the control parameter, here the effective line tension, is increased, the two 
solutions merge and a change of topology is expected at this point. We do not 
know if this change is irreversible since fission requires microscopic reorga- 
nization such as hemifission ( Kozlovsky and Kozlov, 2003| ). Experimentally, 
the daughter vesicles can remain connected by a small filament of lipids but 
if the two vesicles move away, the process is of course not reversible. The 
critical value r c is determined by the fourth equation in EqJHl which gives 
the equilibrium of the forces in the radial direction (axis r). The term in 
t = a/P is due to the line tension and its effect is to pinch the membrane. 
The two others terms (in R\ and R\) are related to the pressure force on the 
membrane and are bounded. The critical value r c is the value for the max- 
imal force on the membrane. For higher line tension (or smaller pressure), 
it is no longer possible to compensate for the line tension which splits the 
system into two independent vesicles. 

One important conclusion of this study is the fact that small domains 
are more easily ejected. This can be validated or invalidated experimen- 
tally when a vesicle has several l Q domains of various size. This conclusion 
is opposite to a floppy membrane whose shape is controlled by elasticity 



(Lipowsky and Dimova, 20031. To show this, we have varied the fraction / 
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of the upper domain (label 1) and we have calculated r c using the data of 
the experimental example (see FigEJ). Since the two phases are equivalent 
when elasticity is neglected, the results are the same for / and 1 — /. The 
parameter r c increases with the size of the smallest domain. This result can 
be explained by a simple argument in the limiting case of a flat domain on 
a flat surface. If the radius of the domain is r, the pinching energy (due 
to the line tension) is approximatively or and the resistance energy (due to 
the pressure) is approximatively Pr 3 . The balance of the two energies gives 
o/P ~ r 2 . Therefore, it is harder to destabilize a large domain than a small 
one. Next, we study the robustness of the model when elasticity is taken into 
account. 



3.2 The elasticity localization. 

Comparing the bending energy (EqJTJ to the osmotic pressure energy (EqJ2J) 
one finds that elastic effect can be neglected if << PR 3 in each phase. 
However, a discontinuity of the tangent appears at the interface between the 
two domains creating a singularity in the curvature. As soon as the bending 
modulus is exactly zero, this discontinuity produces an infinite elastic energy 
contribution, localized near the junction, in contradiction with the weakness 
of elasticity. We are faced with a classical boundary layer model, as found for 
example in the crumpling of an elastic plate JBen Amar and Pomeau,~1 997) 



or the folding of an elastic shell (Pogorelov, 1988). For small but not zero k 



values, near the junction, the elastic effects smooth out the discontinuity by 
locally modifying the shape of the membrane (see Fig. |2|) on a characteristic 



13 



distance of order the elastic length l e in each phase: 

L ■ 

PR, 



Using typical values for giant vesicles (Baumgar t et al., 2003 ), we get l e ~ 
0.5/im, which is very small compared to Ri ~ 10/im. Therefore, we can 
model our system as two spherical caps slightly distorted at the junction on 
a distance of order l e . 

3.2.1 Fold description 

Far away from the fold, the spherical solution (denoted by S) is a good 
approximation but not in the close vicinity of the fold better described by a 
boundary layer (denoted by B) of size l e = l e /Ri- We define a new arclength 
parameter I = (s — §j)/l e and we decompose ■0, r and 7 into 

if) = il*s + iI>b(1) ; r = r s + l e r B (i) ; 7 = 75 + ^7b(0- ( n ) 

with ips = 61 or ips = 62- The quantities ipB, and 75 must vanish far 
away from the junction. Neglecting absorbed molecules, the leading order of 
Eq03 gives 

4> B = simps- (12) 
This is the pendulum equation with solution: 

tan(VW 4 ) = tan(^ cusp /4) exp (±1). (13) 

The plus or minus sign is required for I values, negative or positive: after the 
junction (5 > Sj, I > 0), or before the junction (s < Sj, I < 0). 
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From EqlH and IT2| we derive the elastic energy in each phase: 

Fb ~ - - f ~ I - ibmsfn — Oi 



2nPRf 
+ sin 9i 



L sin 9j < 2 1 — cos 



(14) 



sin 9; — sin 



Ipcusp ~)~ 9j 



The elastic energy (EqiHJ, localized at the junction is proportional to the 
interface length (2iirj = 27rsin6' i in dimensionless parameters) and has the 
same effect as a line tension. Adding the two contributions, we obtain in 
physical units: 



' cusp 



VkiPRi { 2 



cos 



(15) 



+ sin 9i 



\/k 2 PR 2 



+ sin 9 2 



sin 9\ — sin 



cusp 



9i 



sin 9o 



cos • 



sin 



^Pcusp — $2 



+ 



Ipcusp + #2 



The value of ip CU sp is fixed by the boundary conditions (EqEJ 



cusp 



2 arccos 



/R^cos (f ) + v^R^cos (f ) 



(16) 



R\K\ + # 2 k 2 + 2^R 1 R 2 k 1 k 2 cos (^a) 

As expected ^ cusp depends on the ratio of both rigidities. However, it can 
not be easily measured since the size of the fold is very small compared to 
the vesicle size. 

The parameter a cusp measures the strength of elasticity on our spherical- 
cap system. Note that its contribution is angular dependent. Elasticity 
contributes to the line tension and gives an effective line tension a = a + 



a, 



cusp- 



However, the total line tension is now a function of all the physical 



constants (er, P, k q and Kd) which makes it difficult to estimate. Typical 
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values of the elastic line tensions are a cusp ~ 10~ 14 N/m (see Fig. EJ), for 
a ~ lCT 13 iV/m. However, the effect of a cusp on the membrane stability is 
given by the dimensionless number 

"cusp 

n = 

(Pr c - a) 

which measure the relative effect of the elastic contribution with respect to 
the distance at the bifurcation point. In our case, we get an important effect 
with n = 38% and the contribution of the elasticity to the total energy (about 
4%) is not enough to affect the zero-order solution, but can be important for 
the fission of the vesicle. 

3.2.2 Effect on the membrane shape. 

The elastic terms can be taken into account by defining an effective line 
tension. Therefore, the previous results and the capillary solution are still 
valid but with a new control parameter given by f = (a + (? cusp )/P. Note 
that the critical value r c at the bifurcation is still the same. 

A variation of the control parameter f modifies the angles 9\, 9 2 and ip cusp 
and then the elastic line tension a cusp . The figure shows the values of the 
reduced line tension of the fold (a cusp = a cusp /P) versus the reduced total 
line tension f . The solid line is a cusp for the low energy solution. The dashed 
line is a cusp for the high energy solution. The line tension of the fold must 
be smaller than the total line tension a + a cusp since the line tension a due 
to the interface is positive. Therefore, some shapes are no longer physically 
allowed for the unstable solution. The figure [HI shows the energies of the 
vesicle versus the control parameter f for the allowed solutions. 

We have investigated the effect of the size / of the domain 1 on the 
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elastic contribution. In the capillary model, the two domains are equivalent 
and r c is the same for / and 1 — /. The elasticity breaks this symmetry 
since the two domains are no more equivalent: the l domain (here label 2) is 
harder to bend than the la one (here label 1). The figure [HI shows the elastic 
line tension versus the fraction / of the domain 1. Notice that the elastic 
line tension is negative for large domains, meaning that the elasticity fights 
against pinching. 

3.3 Budding by molecule insertion 

The two-cap model remains a solution of the Euler-Lagrange equations when 
molecules are added uniformly. EqlHl connect the concentration of molecules 
to the chemical potential (/Xj) and modify the area Lagrange multiplier 



IH = 2^ - oti(<f>i - (j) eqi ) (17a) 
Hi 



2E<i^ = PR*-2A i <p i R i + a{<f) 2 -<p eq 2 )R 2 (17b) 
PR 2 

ji = — — — sin if> cos if; (17c) 



In a previous paper, we have shown that the two-cap solutions may be unsta- 
ble either above a critical homogeneous concentration given by Ci = PR 2 /Ai 
or for very strong coupling A^/Kj«j >> 1 (All ain and Ben Amar, 2004 ). This 



instability characterizes each phase individually and not the junction between 
phases. Here, we focus on the junction and the experimental conditions are 
assumed to be below these instability thresholds. 
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3.3.1 Fold description 

The interface is the place where strong gradients of molecule distribution are 
found with typical lengthscale given by 

7 - [E 

V OLi 

which must be compared to the vesicle lengthscale Ri. We focus on the case 
where l c « so that concentration gradients are also localized at the fold 
in the elastic boundary layer. For distances larger than l c , the concentration 
of molecules is constant and reaches the value (pi that we choose as unit in 
each phase: so (pi = (pi/ (p. 

Far away from the fold, the sphere (denoted by S) is solution but not in 
the vicinity of the fold, better described by a boundary layer (denoted by B). 
As previously fEq lTT| . we define: 

I = (S ~~ 5 V = ^5+^(0 ; r = r s + ! e r B (l) ; 
7 = 75 + Ljb{1) and = 1 + (p B - 

To describe the fold, we need three dimensionless parameters 



lc ~' l ^mv Xe ~pRU e etXc -^R-j e - (18) 

The conditions for the stability of both phases are A e < 1 and A e A c < 1. 
Expanding the shape equations (EqEJ to leading order gives: 

i>" B = sin ip B + X e (p' B , (19a) 
l 2 c 4>" B = 4>b-\^'b- (19b) 
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The fold energy in the phase i is given by the leading orders of EqSJ 



Fi = 7rl e rjPRf 



ip' B 2 dl — sin ^5 J (sin ip — sin ip s )dl — 2\ e J ip' B dl 



B B B 



+4A e y 4> B dl-2X e J ^ B ^' B dl + j- J 4> 2 B dl + ! 2 c j- J (f)' B 2 dl 

B B C B B 



.(20) 



The energy is the same for both phases. The sum of the two energies is pro- 
portional to rj, the interface length, and defines a new effective line tension 

^ "cusp • 

The uniform insertion of molecules in the vesicle modifies only the La- 
grange multipliers which have no direct physical content, despite the modifi- 
cation of the energy level of the system. The two-spherical cap zeroth order 
solution remains valid without modification of the geometrical parameters 
such as radii and angles at the junction. Therefore, we conclude that the 
bifurcation diagram remains unchanged, except for the values of the energy, 
with the same threshold value found previously. For r < r c , two ideal solu- 
tions still exist, the stable one being observed experimentally. Only gradients 
which appear at the interface modify the cusp shapes and we need to evaluate 
if they are responsible for a change in the line tension value. 

The equations ffTTTf) have no explicit solution but some interesting limits 
can be considered. We focus here on three independent limits: A e << 1, 
A c << 1 and l c « l e . 

First case (A e << 1): the elastic coupling length is small. This limit de- 
couples Eq ll9a,l at zero order, giving exactly the same solution as the case 



without molecule. Eq llQbl allows to calculate the molecule distribution but 
as the terms in A e can be neglected in the energy (Eq l20|) . the effect of the 
molecules is negligible. The elastic line tension is not modified by molecule 
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addition. 

Second case (A c << 1): the weak chemical coupling length. This limit de- 
couples Eq llQbj leading to <pB = <Pb exp (±/// c ) with (f>B the molecule excess 
at the interface, given by the boundary conditions (EqEJ). In physical units 
and taking into account both sides of the fold, we get for the increase of the 
line tension: 



5a cusp = 777 7—0- - (^2 ~ 4>i) 21 

The molecule absorption increases the effective line tension, which may in- 
duce the fission. This effect is only due to chemical gradients near the inter- 
face. It increases with the number of molecules added to the system. 
Third case (l c « 1): the case of a small chemical length. The molecule 
concentration has two very different lengthscales: l c and l e . The chemical 
length l c contributes to the junction between the two domains and can be 
treated as a boundary layer. However, the associated energy is proportional 
to l c and is then negligible. For size larger than l c , EqfTQl becomes: 

4>b = Xc^'b (22a) 
(l-A c A e )V>B = sinV> B (22b) 

and the effective line tension in the phase i is: 



O "cusp 



1 — A e A c ) J ip' B 2 dl — sin ip s J (sin^ — sin ips)dl 



+4A e A c J ip' B dl-2\ e J ip' B dl 



(23) 



B B 

Taking into account the first and second terms leads to the elastic line ten- 



sion (EqUHJ with a multiplicating factor y 1 — A e A c in both phases. Note 
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that A e A c > 1 is not possible in our framework, since the cost associated with 
the concentration gradients (3 are neglected (Allain and Ben Amar, 2004). 



This first contribution indicates that addition of molecules decreases the 
line tension associated with the elastic fold but since it does not depend on 
the concentration it just implies a renormalization of the bending rigidity 



( Leibler, 1986| ). The third term, proportional to A e A c also decreases the line 
tension but does not depend on the molecules concentration. Using physical 
units, the last integral contributes to the line tension by a term proportional 
to the concentration of added molecules: 

Ai0i(0i - Vw) + A 202(Vw - 2 ). 

It is positive when the molecules are inserted in the outer monolayer of the 
membrane (positive A), which is the case found in the experiments. However, 
if the molecules are added in the inner monolayer, it becomes negative and 
budding and pinching are inhibited. In conclusion, the net effect of molecule 
insertion is a decrease of the line tension at low concentration, then a possibly 
increase as the concentration increases if the molecules are inserted from the 
outer monolayer. 

3.3.2 Budding process 

The absorption of molecules does not change the zeroth order shape equations 
of the stretched vesicle. It modifies the shape of the fold near the interface 
giving a new contribution to the effective line tension. If the absorption takes 
place in the external leaflet, it contributes to an increase of the line tension. 
This increase puts the system closer to the bifurcation point controlled by the 
parameter f = (<J+a cusp )/ P and induces a budding of the smaller phase: as f 
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increases, the neck radius decreases (see Fig.GJ)) and the small domain seems 
to lift up. If the concentration is high enough so that f > r c , the budding is 
automatically followed by a fission process, creating two separated vesicles, 
one for each phase. If the concentration is not high enough, the lift-up will 
stop before the change of topology. In the meantime, it is possible that the 
vesicle looses some of its molecules and retracts to its initial configuration. 
This reversibility is impossible when fission is complete for two reasons: first, 
the system relaxes the Gaussian elastic energy and two daughter vesicles may 
be energetically favored, second due to thermal fluctuations, the vesicles move 
away from each other and the coalescence process is unlikely. The fact that 
the fission occurs proves that the time scale for fission is much smaller than 
the possible rearrangement of molecules between the leaflets 

FigureHUreproduces experimental results from Staneva et al. ( |Staneva et al., 2004 1 , 
showing fission of a liquid-ordered domain induced by Phospholipase A 2 pro- 
teins addition. The vesicle is obtained by electro-formation (the electrode 
is visible on the left of the pictures). It is composed by a 45:45:10 mol/mol 
mixture of phosphatidylcholine (PC), sphingomyelin (SM) and cholesterol 
(Choi). The vesicles includes one liquid-ordered domain visible in fluores- 
cence microscopy (not reproduced here): a small fraction (10%) of the PC 
is replaced by a fluorescent lipid analog and is excluded from the l Q phase, 
which appears as a dark circle. The proteins are injected in the neighbor- 
hood of the vesicle by a micropipette (visible on the right of the first picture). 
Phospholipase A 2 activity transforms the PC lipids into LysoPC, a conical 
molecule, by cutting one of the two hydrophobic tails. Fission occurs about 
10 seconds after protein injection. 
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Similar fission process have been observed when detergents like LysoPC, 
Triton XI 00 or Brij 98 are added in important quantities near a similar vesicle 



(St aneva et al., submitted! ) . However, in this case, the fission is not always 
complete: the daughter vesicles may remain connected by a small lipid fila- 



ment, as also observed in (Tanaka et al., 2004 1. This is not in contradiction 
with our model since the fission process requires to split the lipid bilayers 
at the molecular level, which is out of reach of our treatment. This level 
requires a microscopic description as done in ( |Kozlovsky and Kozlov, 20 03). 
The fission process prefers small domains, as predicted by our model. If the 
concentration in detergents is not high enough or if the Phospholipase A 2 is 
not activated, the liquid-ordered domains bud without complete fission. It is 
also possible to observe a relaxation of the vesicle, which recovers its initial 
shape. 



4 Conclusion 

Our model explains why ejection of a domain from an inhomogeneous vesicle 
can be achieved by osmotic shocks or molecule absorption. It is based on 
physical stability concepts in the spirit of the existence and stability anal- 
ysis of the well-known catenoid. We predict a complete irreversible fission 
above a critical parameter. From a macroscopic point of view, the complete 
fission is favored, it decreases the total energy of the system at the threshold 
of stability because of the Gaussian energy. This fission can be inhibited if 
a membrane thread exists between the two phases. The existence of such 
a thread is out of reach of our approach. If it does not exist, the vesicles 
separate from each other. If it exists and if the experimental forcing relaxes, 
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the two vesicles may fuse in principle. The experiments discussed here are 
in favor of a complete fission mechanism. For simplicity, the model is re- 
stricted to two domains of different sizes: extension to multi-phase domains 
complicates the geometry but will not change the physical results. 

Acknowledgments 

We would like to thank M. Angelova, P. Bassereau, T. Baumgart and G. 
Staneva for various discussions on experimental aspects and A. Goriely for 
critical reading. 

References 

Allain, J.-M. and M. Ben Amar. 2004. Biphasic vesicle: instability induced 
by adsorption of proteins. Physica A. 337:531-545. 

Allain, J.-M., C. Storm, A. Roux, M. Ben Amar and J.-F. Joanny. 2004. 
Fission of a multiphase membrane tube. Phys. Rev. Lett. 93:158104. 

Baumgart, T., S. T. Hess and W. W. Webb. 2003. Imaging coexisting fluid 
domains in biomembrane models coupling curvature and line tension. 
Nature. 425:821-824. 

Ben Amar, M., P. P. da Silva, N. Limodin, A. Langlois, M. Brazovskaia, 
C. Even, I. V. Chikina and P. Pieranski. 1998. Stability and vibrations 
of catenoid-shaped smectic films. Euro. Phys. J. B. 3:197-202. 

Ben Amar, M. and Y. Pomeau. 1997. Crumpled paper. P ROY SOC LOND 
A MAT 453:729-755. 

24 



Bickel, T., C. Jeppesen and C. M. Marques. 2001. Local entropic effects of 
polymers grafted to soft interfaces. Eur. Phys. J. E. 4:33-43. 

Brown, D. A. and E. London. 2000. Structure and function of sphingolipid- 
and cholesterol-rich membrane rafts. J. Biol. Chem. 275:17221-17224. 

Chazal, N. and D. Gerlier. 2003. Virus entry, assembly, budding, and mem- 
brane rafts. Microbiol. Mol. Biol. Rev. 67:226-237. 

de Almeida, R. F. M., A. Fedorov and M. Prieto. 2003. Sphin- 
gomyelin/phosphatidylcholine/cholesterol phase diagram: Boundaries 
and composition of lipid rafts. Biophys. J. 85:2406-2416. 

Dietrich, C, L. A. Bagatolli, Z. N. Volovyk, N. L. Thompson, M. Levi, 
K. Jacobson and E. Gratton. 2001. Lipid rafts reconstituted in model 
membranes. Biophys. J. 80:1417-1428. 

Dobereiner, H.-G., E. Evans, M. Kraus, U. Seifert and M. Wortis. 1997. Map- 
ping vesicle shapes into the phase diagram: A comparison of experiment 
and theory. Phys. Rev. E. 55:4458-4474. 

Helfrich, W. 1973. Elastic properties of lipid bilayers: Theory and possible 
experiments. Z. Naturforsch A. 33:305-315. 

Holopainen, J. M., A. J. Metso, J.-P. Mattila, A. Jutila and P. K. J. Kin- 
nunen. 2004. Evidence for the lack of a specific interaction between 
cholesterol and sphingomyelin. Biophys. J. 86:1510-1520. 

Jaric, M., U. Seifert, W. Wintz and M. Wortis. 1995. Vesicular instabilities: 
The prolate-to-oblate transition and other shape instabilities of fluid 
bilayer membranes. Phys. Rev. E. 52:6623-6634. 

25 



Jtilicher, F. and R. Lipowsky. 1996. Shape transformations of vesicles with 
intramembrane domains. Phys. Rev. E. 53:2670-2683. 

Kohyama, T., D. M. Kroll and G. Gompper. 2003. Budding of crystalline 
domains in fluid membranes. Phys. Rev. E. 68:061905. 

Kozlovsky, Y. and M. M. Kozlov. 2003. Membrane fission: Model for inter- 
mediate structures. Biophys. J. 85:85-96. 

Laradji, M. and P. B. Sunil Kumar. 2004. Dynamics of domain growth in 
self-assembled fluid vesicles. Phys. Rev. Lett. 93:198105. 

Leibler, S. 1986. Curvature instability in membranes. J. Phys. (France). 
47:507-516. 

Li, X.-M., M. M. Momsen, J. M. Smaby, H. L. Brockman and R. E. Brown. 
2001. Cholesterol decreases the interfacial elasticity and detergent solu- 
bility of sphingomyelins. Biochemistry. 40:5954-5963. 

Lipowsky, R. and R. Dimova. 2003. Domains in membranes and vesicles. J. 
Phys.: Condens. Matter. 15:S31-S45. 

Lipowsky, R. and E. Sackmann. 1995. Structure and Dynamics of Mem- 
branes, Handbook of Biological Physics. Elsevier, North-Holland. 

Pogorelov, A. 1988. Bending of surfaces and stability of shells. Volume 72, 
Translations of mathematical monographs. 

Seifert, U. 1993. Curvature-induced lateral phase segregation in two- 
component vesicles. Phys. Rev. Lett. 70:1335-1338. 



26 



Seifert, U. 1997. Configurations of fluid membranes and vesicles. Adv. in 
Phys. 46:13-137. 

Sens, P. 2004. Dynamics of nonequilibrium membrane bud formation. Phys. 
Rev. Lett. 93:108103. 

Siegel, D. P. and M. M. Kozlov. 2004. The gaussian curvature elastic modulus 
of n-monomethylated dioleoylphosphatidylethanolamine: Relevance to 
membrane fusion and lipid phase behavior. Biophys. J. 87:366-374. 

Simons, K. and E. Ikonen. 1997. Functional rafts in cell membranes, nature. 
387:569-572. 

Staneva, G., M. I. Angelova and K. Koumanov. 2004. Phospholipase a(2) 
promotes raft budding and fission from giant liposomes. Chem. Phys. 
Lipids 129:53-62. 

Tanaka, T., R. Sano, Y. Yamashita and M. Yamazaki. 2004. Shape changes 
and vesicle fission of giant unilamellar vesicles of liquid-ordered phase 
membrane induced by lysophosphatidylcholine. Langmuir. 20:9526- 
9534. 

van Meer, G. and H. Sprong. 2004. Membrane lipids and vesicular traffic. 
Curr. Opin. Cell Biol. 16:373-378. 

Veatch, S. L. and S. L. Keller. 2003. Separation of liquid phases in giant 
vesicles of ternary mixtures of phospholipids and cholesterol. Biophys. 
J. 85:3074-3083. 



27 



28 



A 



4> 



r(s) 



Figure 1: General parameterization of an axisymmetric vesicle in cylindrical 
coordinates. The dashed curve is the membrane. The parameterization is 
done by the arc-length s. The shape of the membrane is given by r(s) and 
ip(s). The two domains have the same parameterization. 
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Figure 2: Schematic representation of a axisymmetric vesicle, including the 
four parameters Ri, R 2 , Q\ and 9 2 used in the vesicle description. The circle 
details the fold near the interface, where the elastic properties can no longer 
be neglected. 
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Figure 3: Black and white version of figure (lb) from Baumgart et al.'s 



work ( |Baumgart et al., 2003| ). The picture is a two-photon microscopy im- 
age, showing equatorial section of GUVs with two coexisting domains. The 
l Q domain appears in grey here and the Id in dark. Scale bar, 5 \im. Repro- 



duced from (Baumga rt et al., 2003 ) with the authorizations of the authors 
and editor. 
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a) b) 




Figure 4: The two solutions of the E.-L. equations (EqEI and EH) for r = 
20.5/im 2 . The two pictures have the same scale, (a) Experimental solution 



( Baumgart et al., 2003D , with Ri = 5.3fim, R 2 = llfxm, Q\ = 1.3 and 9 2 = 
0.51. The associated dimensionless energy, given by EqUHlis Ftot = —1.380. 
(b) Calculated solution with Ri = 4.0/zm, R 2 = 10/im, 61 = 2.0 and 9 2 = 
0.36. The dimensionless energy of vesicle (b) is Ftot = —1.377, meaning 
that the solution is experimentally unstable. 
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Figure 5: Dimensionless energies Ftot (figure a) and interface radius rj 
(figure b) of the solutions of the E.L. equations (EqEI and versus the 
control parameter r = a /P. The calculation has been done with the area 
A\ = 136/xm 2 and A 2 = 1296/im 2 . In both figures, the solid line corresponds 
to the stable solution, experimentally observable, the dashed line to the un- 
stable solution and the dotted line to f = f c , the critical value of the control 
parameter. For f > f c , there is no longer a solution. Four pictures of vesicles 
showing the shape transformation with r have been added. 
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Figure 6: Values of the critical control parameter r c versus the fraction / 
of the upper domain (label 1). The calculation has been done with a fixed 
total area A to t = A\+ A 2 = 1433/im 2 . The areas of the domains are given by 
A\ = f A to t and A% — (1 — f)A to t- The dashed line shows the fraction 0.5. 
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Figure 7: Reduced fold line tension a cusp versus the reduced total line tension 
f for the two possible solutions of EqlHl The solid line corresponds to the 
stable solution, the dashed line to the unstable solution. The fold line tension 
has been calculated using EqlT5l and fT6l The parameters are A\ = 136/xm 2 , 
A 2 = 1296/im 2 (the fixed area of each phase), P = lO~ 2 Pa, K\ = 10~ 19 J 
and k 2 = 1CT 18 J. The dotted line separates the possible solutions from the 
impossible one. For a cusp above this line, the line tension associated to the 
fold is greater than the total line tension, requiring a negative line tension at 
the junction between the two domains, which is impossible. 
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Figure 8: Dimensionless energies F TO t of the solutions of the E.L. equa- 
tions (EqEI and IHJ) versus the control parameter f including the effect of the 
elastic fold. The parameters are A\ = 136/im 2 , A 2 = 1296/xm 2 , P = lCT 2 Pa, 
Ki = 10 _19 J and n 2 = 10 _18 J. The solid line corresponds to the stable 
solution, the dashed line to the unstable solution, the dotted line to f = r c , 
the critical value of the control parameter. For f > r c , there is no longer a 
solution. 
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Figure 9: Fold line tension a cusp for the critical value of the total line tension 
f = f c , versus the fraction / of the liquid-ordered domain. The calculation 
has been done with a fixed total area A tot = A\ + A 2 = 1433/im 2 , the other 
parameters being P = lCT 2 Pa, K\ = 10~ 19 J and n 2 = 10 _18 J. The areas of 
the domains are given by A\ = f A tot and A 2 — (1 — f)A tot . The dashed line 
shows the fraction 0.5 and the dotted line a cusp = 0. The l Q and Id domains 
(resp. label 2 and 1) do not have the same effect since the elastic moduli are 
not equal. Contrary to a true line tension, this effective line tension can be 
negative for / near 0.5. 
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Figure 10: Angles 9i, 9 2 and ip cusp for f = f c , the critical value of the 
control parameter, versus the fraction / of the liquid-ordered domain. The 
calculation has been done with a fixed total area A to t = A\ + A 2 = 1433/xm 2 , 
the other parameters being P = lCT 2 Pa, K\ = 10 _19 J and k 2 = 1CT 18 J. 
Label 1 corresponds to the Id domain and label 2 to the l a . The areas of the 
domains are given by A\ = f A to t and A 2 = (1 — f)A to t- The stars are for 
the angle 9%, the circles for 9 2 and the triangles for ip C us P - The angle ipcusp 
is always closer from the angle 9 2 since the liquid-ordered domain is hard to 
bend. 
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Figure 11: Ejection of a liquid-ordered domain induced by proteins Phospho- 
lipase A2. The domain is visible on fluorescence microscopy (not reproduced 
here). The proteins are injected by micropipette (figure a). The liquid- 
ordered domain buds (figure b and c) before the fission (figure d). Repro- 



duced form flStaneva et al., 2004D with the authorization of the editor. Bar: 
20 fim. 
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